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Abstract 

We present the results of the current analysis of the partial wave 

IjPC = 00++ based 

on the available data for meson spectra (tttt, KK, r]T], T]r]' , mnnr) . In the framework of the 
' ii'-matrix approach, the analytical amplitude has been restored in the mass region 280 

■ MeV< ^/s < 1900 MeV. The following scalar-isoscalar states are seen: comparatively nar- 

^ ; row resonances /o(980), /o(1300), /o(1500), /o(1750) and the broad state /o(1200 - 1600). 

' The positions of the amplitude poles (masses and total widths of the resonances) are deter- 

mined as well as pole residues (partial widths to meson channels tttt, KK ,r]rj,rjrj' jTrrnnr). 
, The fitted amplitude gives us the positions of the iC-matrix poles (bare states) and the 

I values of bare-state couplings to meson channels thus allowing the quark-antiquark nonet 

• classification of bare states. On the basis of the obtained partial widths to the channels 

TTTT, KK ,ri7],r]r]' , we estimate the quark/gluonium content of /o(980), /o(1300), /o(1500), 
PhI /o(1750), /o(1200 - 1600). For /o(980), /o(1300), /o(1500) and /o(1750), their partial 

widths testify the qq origin of these mesons though being unable to provide precise eval- 
uation of the possible admixture of the gluonium component in these resonances. The 
ratios of the decay coupling constants for the /o(1200 — 1600) support the idea about 
^ ' gluonium nature of this broad state. 



1 Introduction 



The classification of meson states in the scalar-isoscalar sector is a key problem for understand- 
ing the strong QCD being a subject of intensive discussions in recent years, see, for example, 
|l|, 1^, 1^, ^, 1^, 1^, 1^ and references therein. 

In this paper we present the results of the ii'-matrix analysis of the 00++ wave in the invariant 
mass range 280-1900 MeV. This analysis is a continuation of earlier work ^, |I0|, |1T|. In 

the 00++ wave had been reconstructed on the basis of the following data 

1^1) GAMS data on the 5'-wave two-meson production in the reactions Tcp —>■ n^n^n, rjrjn and 



the latter paper []11 
set: 
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7]r]'n at small nucleon momenta transferred, \t\ < 0.2 (GeV/c)^ ||T2| , 13 



(2) GAMS data on the vrvr S-wave production in the reaction np TT^ir^n at large momenta 
transferred, 0.30 < |t| < l.O jGeV/c)^ p[; 

(3) BNL data on np- KKn Q; 

(4) CERN-Miinich data on tt+tt- tt+tt- ^ 



(5) Crystal Barrel data on pp{dX rest, from liquid H2)—^ tt^tt^tt^, ir^TT^r], ir^r]'!] |16 . 

Now the experimental basis has been much broadened, and additional samples of data are 
included into present analysis of the 00++ wave, as follows: 

(6) Crystal Barrel data on proton-antiproton annihilation in gas: pp{a.t rest, from gaseous 

(7) Crystal Barrel data on proton-antiproton annihilation in liquid: pp{at rest, from liquid 
H2)^ TT+TT-vrO, K+K-TT^, KsKsTT^, K+KsTT- p, 0; 



(8) Crystal Barrel data on neutron-antiproton annihilation in the liquid deuterium: r;,p(at rest, 
from liquid Da)^ 7r°7r°7r-, 11-^-71+, KsR-tt^, KsKg-n" [Q, 0; 

(9) E852 Collaboration data on the tttt S-wave production in the reaction vr^p 7r°7r°n at the 
nucleon momentum transfers squared < |t| < 1.5 (GeV/c)^ ]TP[ . 



The production of resonances in the 00++ wave is accompanied by a considerable back- 
ground. This circumstance forces us to carry out a simultaneous analysis of as much as possible 
number of spectra, for only in such an analysis one can single out resonance parameters reliably. 



As compared to the paper |TT|, the reactions of the pp annihilation in gas have been included 
into the present analysis. One should keep in mind that, while in the liquid hydrogene the pp 
annihilation is going dominantly from the S-state, in the gas there is a large admixture of the 
P-wave, thus providing us an opportunity to analyse three-meson Dalitz plots in more detail. 

New Crystal Barrel data allow us to study the two-kaon channel more reliably as compared 
to what had been done before. This is undoubtedly important for the conclusion about the 
quark-gluonium content of the /o-mesons under investigation. 

The data of the E852 Collaboration on the reaction -K^p -K^-K^n at piab = 18 GeV/c 



19[] together with GAMS data on the rection vr^p n^-n^n at piab = 38 GeV/c give us 
a good opportunity for the study the resonances /o(980) and /o(1300), for at large momenta 
transferred to the nucleon, \t\ ~ 0.5 ~ 1.5 (GeV/c)^, the production of resonances goes with 
a small background thus allowing us to fix convincingly their masses and widths. This is 
especially important for /o(1300): in the compilation this resonance is referred as /o(1370), 
with the mass in the interval 1200 — 1500 MeV, though experimental data favour the mass near 
1300 MeV. 

In Section 2 we present a set of formulae for the analysed i^'-matrix amplitudes as follows: 

(i) the K-maXiYx. amplitudes of the mass-on-shell transitions into channels tttt, KK, rir], rir]', inrmr] 

(ii) the /T-matrix amplitudes for the transitions irir^^^t) tttt and 7ra[^\t) tttt, where 
7r^^\t) and a[^\t) are reggeized pion and ai-meson in the high energy reaction n^p —>■ (7r7r)s n; 

(iii) the iT-matrix three-meson production amplitudes in the pp and np annihilations. 

The fC-matrix amplitude determines both the amplitude poles (masses and widths of res- 
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onances) and the ii'-matrix poles (masses of bare states). The i^-matrix poles differ from the 
amplitude poles in two points: 

(i) The states corresponding to the i^'-matrix poles do not contain any component associated 
with real mesons which are inherent in resonances. The absence of a cloud of real mesons allows 
us to refer conventionally to these states as bare ones ^ |lT|. 

(ii) Due to the transitions bare state(l) real mesons — > bare state(2) the observed resonances 
are the mixtures of bare states. So, for the quark systematics, bare states are the primary 
objects rather than resonances. 

In Section 3, we present the main results of the present analysis. As in previous papers, 
we have had three solutions which are denoted, following Solutions I, II-l and II-2. In 



all Solutions, the characteristics of the resonances /o(980), /o(1300), /o(1500) and broad state 
/o(1200 — 1600) coincide with each other and coincide, with a reasonable accuracy, with what 
had been obtained in From this point of view, the performed analysis provided us with 
unambigously determined properties of physical states, with an exception for the resonance 
/o(1750), for which the channel /o(1750) — *• Tnnm is not well defined, that resulted in different 
full widths given by Solutions I and II. Solutions I, II-l and II-2 differ by the characteristics 
of bare states and background terms in the iiT-matrix, just as it was the case in Inl, though 



it must be particularly emphasized that the present analysis reveals a noticeable tendency for 
the parametrs from different Solutions to become closer to each other. In the present analysis 
the distinction in characteristics of bare states is not dramatic with an exception, again, for 
/o(1750), for which the coupling constant of the i^- matrix pole to the channel vrTrvrvr given by 
Solutions I and II have opposite sign. 

The relations between coupling constants of the qq bare state, /q°^*^, to channels vrvr, KK, rjr], 
rjrj' provide us with an information on the relative weight of the components nn = {uu + dd) / 
and ss. For the gluoniun, the relations between coupling constants are nearly the same as for 
the flavour qq singlet, thus leading to the ambiguity in the definition of the glueball: Solutions 
I, II-l and II-2 provided us with different possibilities for the glueball among bare states. 

In Solution I, there are two bare states, every of them pretending to be the glueball: 
/o^«''^(1220 ± 30) and /o^'^'^^(1635 ± 25). In Solution II-l only one state, /o^"'-^(1225 ± 25), satis- 
fies the constraint inherent in the glueball. In Solution II-2, again there are two states whose 
coupling constants are appropriate to the glueball: /q"™(1230 ± 30) and /g"'"^(1560 ± 25) . The 
uncertainty in the classification of bare states, which is based on the value of a coupling to 
channels /o"*^*^ t^tt, KK, rjr], rjr]', makes us to apply the other methods to single out the states 
which are extra ones for the gg-systematics, namely, exotic states. The method based on the 
consideration of trajectories in the (n, M^)-plane (n is radial quantum number of the qq state 
and M its mass) is discussed in Section 5. 

On the basis of the extracted partial widths for the channels vrvr, KK, rjr], tjt]', the following 
properties of resonances /o(980), /o(1300), /o(1500), /o(1750) and broad state /o(1200 - 1600) 
are to claim: 

1. /o(980): This resonance is dominantly the qq state, qq = nnsiiup + ss cos (p, with a 
large ss component. Under the assumption that the admixture of the glueball component is 
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not greater than 15%, Wgiuonium ^ 15%, the hadronic decays give us the following constraints: 
—95° < < —40°, that is in agreement with data on radiative decays /o(980) — >• 77 and 
0(1020) — > 7/0(980) [^, Rather large uncertainties in the determination of mixing angle 



are due to the sensitivity of coupling constants to plausible small admixtures of the gluonium. 
If the gluonium component is absent, hadronic decays provide ip = —67° ± 10°. 

As concern the characteristics of the /o(980), we cannot observe any special property which 
would single it out from a qq series: it belongs to linear qq trajectory on the {n, M^) and 
(J, M^) planes and its decay couplings to hadronic channels are of the same order as for the 
other /o resonances, see Table 5 in Section 4. Besides, this resonance is produced, without 
any suppression, in hard processes as well as in heavy-quark decays, e.g. in the reactions 
n^p (vrTr)^ n at large |t| |12, 19| and tt+tt^tt^ [23|. 



2. /o(1300): This resonance is the descendant of the bare qq state which is close to the 
flavour singlet. The resonance /o(1300) is formed due to a strong mixing with the primary 
gluonium and neighbouring qq states. The qq content, qq = nncos(p + sssimp, in the /o(1300) 
strongly depends on the admixture of the gluonium component. At Wgiuonium ^ 30% the 
mixing angle changes, depending on Wgiuonium, in the interval —45° < v5[/o(1300)] < 25°; at 
Wgiuonium = the hadronic decays provide </9[/o(1300)] = —7° ± 10°. 

3. /o(1500): The resonance is the descendant of the bare state with a large nn component. 
Being similar to /o(1300), the resonance /o(1500) is formed by mixing with the gluonium as 
well as with neighbouring qq states. The qq content, qq = nn cos (p + ss sin p, depends on 
the admixture of the gluonium: at Wgiuonium ^ 30% the mixing angle changes, depending on 
Wgiuonium, in the interval -20° < ip[fo{1500)] < 25°; at Wgiuonium = one has Lp[fo{1500)] = 
10° ± 10°. 

4. /o(1750): This resonance is the descendant of the bare state of the radial excitation 
nonet 2^Piqq, which flavour wave function has large ss component. In Solutions I and II, the 
resonance /o(1750) has different values of the ss component. In Solution I, the ss component 
dominates: in the absence of gluonium component, <y9[/o(1750)] = — 72°±5°, and if the gluonium 
admixture is Wgiuonium ^ 30%, then —110° < (p[/o(1750)] < —35°. In Solution II, in the 
absence of gluonium component, v9[/o(1750)] = —27° ± 10°, and with gluonium admixture 
Wgiuonium < 30%, -50° < <^[/o(1750)] < 10°. 

5. /o(1200 — 1600): The broad state is the descendant of a primary glueball. The analysis of 
hadronic decays of this resonance confirms its glueball nature, for the qq component is allowed 
to be in a state produced by the glueball: {qq)giuebau = {\/2nn + y/X ss)/y/2 + A where 



A defines relative probability to produce the ss pair by gluon field, uu : dd : ss = 1 : 1 : \. 
Different estimations give A ~ 0.5 - 0.8 % ^ thus leading to (^[/o(1200 - 1600)] ~ 27° - 32°. 
Just such a magnitude appears in hadronic decays of the broad state, though the value of a 
possible admixture of the {qq)giuebau cannot be fixed by hadronic decays. The impossibility to 
determine the quark-antiquark component is due to the fact that the relations between decay 
coupling constants are exactly the same for the gluonium and {qq)giuebaii- 

Determination of the /o(1200 — 1600) and its characteristics needs additional comments. 
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The spectra in the channels tttt, KK, rjr], rjr]' require the introduction of a broad bump, and 
this bump occures to be universal in different reactions thus making it possible to describe 
it as a broad resonance. One of a distinct characteristics of resonances is a factorization of 
the resonance amplitudes: the amplitude, gin{s — M'^)~^gout, may be factorized with universal 
couplings {gin, gout) for a variety of reactions. The fitting to a large number of reactions 
studied here gets along with the factorization of the broad state. Large width of the braod 
state does not allow us to fix reliably the mass of /o(1200 — 1600), but its strong production in 
a variety of reactions under investigation makes it possible to assuredly determine the ratios of 
couplings to the channels tttt, KK, rjr], rjr]' thus determining reliably the gg/gluonium content 
of /o (1200 - 1600). 

The present i^-matrix analysis does not point to the existence of a light a-meson which is 
actively discussed now (see e.g. ||l]] and references therein), in particular in connection with the 
recently reported signals in tt+tt+tt" with the amplitude pole at M = (480±40) —i{160± 

30) MeV 0, in J/^ nnuj with the pole at M = (390 tfe) - ^(141 1|) MeV and in 
r -R-R-KV with the pole at M ~ 555 — ^270 MeV [^. Possible explanation of this contradiction 
may lay in a strong suppression of the light cr-meson production in the annihilation process 
pp — ^ TTTTTT, though there is no visible reason for this suppression (recall that the statistics for the 
Crystal Barrel reactions is by two orders of magnitude larger than in P3| , P^ , P7[|). Alternative 
explanation can be related to a restricted validity of the i^-matrix approach at small invariant 
energy squared s: due to the left-hand cuts in the partial amplitude, which were not properly 
taken into account in the i^-matrix approach, one may believe that the analysis does not 
reconstruct anlytical amplitude in due course at Res < 4m^, i.e. at (ReM) — (ImM)^ < 4m^. 
In a number of analyses, including those performed in the dispersion relation technique where 
the left-hand cut can be accounted for in one way or another, the pole ascribed to the a-meson 



had been obtained just at Re s < 4m^, see e.g. ^ ^ ^ 



In Section 3, to illustrate the level of accuracy for the Dalitz-plot description in the reactions 
pp three mesons and np three mesons we demonstrate experimental spectra in various 
reactions versus fitted curves of Solution II-2. We also show the tttt spectra in the reactions 
TT~~p {TnT)sn measured by GAMS [0 and E852 |jl9[ and their description in one of the fits 



(a detailed discussion of these spectra is presented in a separate article |^). One important 
statement of this analysis is worth being stressed here: for a combined description of the GAMS 
and E852 data, one needs to introduce the t-channel exchanges by the leading and daughter 

Regge trajectories: irieading, 0'l{leading) and T^daughter, 0'l(daughter)- 

The parameters of the i^-matrix amplitude obtained directly from the fit (they are given 
in Tables 2, 3, 4 for Solutions I, II-l and II-2) serve us as characteristics of bare states. To 
define coupling constants of the /o-resonances to channels tttt, KK, rjr], rjrj', utititi one needs 
to calculate the amplitude residues, see Section 4. In this Section, on the basis of calculated 
couplings, we give partial widths of resonsnces /o(980), /o(1300), /o(1500) and /o(1750). 

In Section 5, the calculated couplings are analysed in terms of quark combinatorics rules to 
define the quark-antiquark and gluonium content of resonances. 

The analysis performed here does not give us a precise correlation between the nn, ss and 
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gluonium components. The reason of this is not a possible insufficient accuracy of data but the 
structure of /o-mesons themselves in the presence of gluonium admixtures — this problem is 
emphasized in Conclusion. 

In Appendices A, B, C we present, in more detail, the formulae used in this analysis. 



2 [IJ = 00++)-wave: the i^^-matrix amplitude and quark- 
combinatorics relations for the decay couplings 

Here we set out the ii'-matrix formulae for the analysis of the (/ J^*" = 00^"'')-wave and present 
the quark-combinatorics relations for the i^-matrix couplings which are used for the study the 
gg/gluonium content of the states under consideration. 



2.1 i^'- matrix amplitude 

The i^T-matrix technique is used for the description of the two-meson coupled channels: 

A = k{i ~ipk)~\ (1) 

where K is n n matrix (here n is the number of channels under consideration) and / is 
unit matrix. The phase space matrix is diagonal: pah = SabPa- The phase space factor pa is 
responsible for the threshold singularities of the amplitude: to keep the amplitude analytical in 
the physical region under consideration we use analytical continuation for pa below threshold. 
For example, the 7777 phase space factor p^^ = (1 — 4m^/s)^/^ is equal to z(4m^/s — 1)^/^ below 
the rji] threshold (s is the two-meson invariant energy squared). To avoid false singularity 
in the physical region, we use for the rjr]' channel the following phase space factor p^j^/ = 
(1 - (m^ + m^/)Vs)i/2. 

For the multi-meson phase volume, we use the four-pion phase space defined as phase space 
of either pp or crcr, where a denotes the S'-wave vrvr amplitude below 1.2 GeV. The result does 
not depend practicaly on whether we use pp or aa state for the description of multi-meson 
channel: below we write formulae and the values of obtained parameters for the pp case, for 
which the fitted expressions are less cumbersome. 



2.1.1 Scattering amplitude 

For the S'-wave scattering amplitude in the scalar-isoscalar sector, we use the parametrization 
similar to that of [|, |10|, |ll|]: 

KAs)=rEj^^+U ^^^^ , (2) 
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where K^j^ is a 5x5 matrix {a,b = 1,2,3,4,5), with the following notations for meson states: 
1 = TTTT, 2 = KK, 3 = rji], 4 = rjr]' and 5 = multimeson states (four-pion state mainly at 
v/i < 1.6 GeV). The is the coupling constant of the bare state a to meson channel; the 
parameters fab and sq describe a smooth part of the ii"-matrix elements (1 < sq < 5 GeV^). 
We use the factor (s — sa)/{s + sao) to suppress false kinematical singularity at s = in the 
physical region near the vrvr threshold. The parameters sa and sao are kept to be of the order 
of Sa ~ (0.1 — 0.5)m^ and s^o ~ (0.1 — 0.5) GeV^; for these intervals, the results do not depend 
practically on precise values of sa and s^o- 

For the two pseudoscalar-particle states, vrvr, KK, rjr], rjrj' , the phase space matrix elements 
are equal to: 

pa{s) = ]J , a = 1,2, 3, 4 (3) 

where mia and are the masses of pseudoscalars. The multi-meson phase-space factor is 
determined as follows: 

. . f Psi at s < 1 GeV^, , 
P^^'^ = [Z.ts>l GeV< 

dsi f ds2 



P51=P0 I -^J ^M^r(si)r(s2) V(S + Si- - 4SS, 



X 



( s - lQml\ 

Here si and S2 are the two-pion energies squared, M is the p-meson mass and T{s) its energy- 
dependent width, r(s) = 7pf (s). The factor po provides the continuity of p^^s) at s = 1 GeV^. 
The power parameter n is taken to be 1, 3, 5 for different variants of the fitting; the results are 
weakly dependent on these values (in our previous analysis |TT| the value n = b was used). 



2.1.2 High-energy production of the S'-wave mesons vrvr, KK, rjr] and rjr]' in the vrp 
collisions 

Here we present the formulae for the high-energy S'-wave production of vrvr, KK, rjr], rjr]' at 



small and moderate momenta transferred to the nucleon. In [|T^, |T3[ 14, 19|, the np collisions 
were studied at pbeam ~ (15 — 40) GeV/c (or s.j,n — '^iriNPbeam ~ 30 — 80 GeV^). At such 
energies, two pseudoscalar mesons are produced due to the t-channel exchange by reggeized 
mesons belonging to the vr and Oi trajectories, leading and daughter ones. 

The vr and ai reggeons have different signatures, C,n = +1 and C,ai = —1- Accordingly, we 
write the vr and ai reggeon propagators as 

sin(|a^(t)) ' cos(faai(i)) 
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Following we use for leading trajectories: 

a^ileading) (t) ^ -0.015 + 0.72t, aameading) [t) ^ "0.10 + 0.72t, (6) 

and for daughter ones: 

a7T{daughter){t) - -1-10 + 0.72t, aal{daughter){t) - "I-IO + 0.72t . (7) 

Here the slope parameters are in (GeV/c)~^. In the centre-of-mass frame, which is the most 
convenient for the consideration of reggeon exchanges, the incoming particles move along the 
z-axis with a large momentum p. In the leading order of the 1/p expansion, the spin factors 
for TT and ai trajectories read: 

TT — trajectory : , (8) 

ai — trajectory : i{(^nz) , 

where = p/p and q± is the momentum transferred to the nucleon {t ~ —(1±)- The Pauli 
matrices a work in the two-component spinor space for the incoming and outgoing nucleons: 
{v*out^'^in) (for more detail see, e.g., @, ^). Consistent removal from the vertices (H) of terms 
decreasing with p — > oo is necessary for the correct account for daughter trajectories which 
should obey, similarly to leading ones, the constraints imposed by the t-channel unitarity. 

In our calculations, we modify reggeon propagators in (^) conventionally by replacing 

s-kn — > — , (y) 

S-kNO 

where the normalization parameter s^^nq is of the order of 2-20 GeV^. To eliminate the poles 
at t < we introduce additional factors (Gamma-functions) into reggeon propagators by sub- 
stituting in (^ as follows: 

sin i'^a^it)] ^ sin f^a^(t)') V { ^ + 1^ , (10) 



2 V V2 

cos {^aalit)^ cos (^|«al(t)) T + 1^ . 

The i^-matrix amplitude for the transitions nR{t) tttt, KK, tjt], rji]', tttttctt, where R{t) refers 
to reggeon, reads: 

A^n = KMi-^pKr\ (11) 
where K^r is the following vector: 

Here G^^^{t) and Fj^ji^blt) are the reggeon t-dependent form factors. The following limits are 
imposed on the form factors: 

[t ml) = gi^J , F..,,(t ^ ml) = , (13) 

where gl^ and fTrn,b enter the matrix element (Q). 
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2.2 Three-meson production amplitudes 



In this Section we present the formulae for the reactions pp — > 7r°7r°7r°, tt^tt^t;, ir'^rjr] from the 
hquid H2, when annihilation occurs from the ^Sopp state and scalar resonances, /o and Qq, 
are formed in the final state only. Of course, this is only one sub-process from a number of 
reactions considered here; still, it is rather representative with respect to the applied technique 
of the three-meson production reaction. A full set of amplitude terms taken into account in the 
present analysis (production of vector and tensor resonances, pp annihilation from the P-wave 
states ^Pi, ^P2, ^Pi) is given in Appendices A and B. 

For the transition pp (^Sq) — > 7r°7r°7r°, the amplitude has the following structure: 



Ipp (iSo)-»7r07r07rO 



^75 



-^(?i) 



(14) 



N 



The amplitude (is'o)7rO,7r07ro('Sij) describes the diagrams with interacting mesons, with the last 
interaction for the particles i and j. In the initial state, the four-spinors ■0(— ^2) and ■0(?i) refer 
to antiproton and proton, respectively. 

The amplitudes for the transitions pp (^Sq) — > rjir^TT^, n^'qr] have similar form: 



and 



^PP (i5o)-»»?7r07rO 



V'(-?2); 



X 



^pp(i5o)r),7r07ro(s23) + ^pp (i5o)7rO,r)7rO (^13) + ^pp (i5o)7rO,r)7rO (^12) 



(15) 



?75 



2\/2mM 



X 



App (i5o)7rO,7?r)(s23) + ^pp (i5o)T),??7rO (^13) + ^pp (i5o)r),7?7rO (^12) 



(16) 



The following amplitude is used for the two-meson interaction block in the scalar-isoscalar state 
{IJ — 00) in the reactions pp tt^tt^tt^ and pp — > Tr^rjr]: 



-^pp (i5o)7rO,6('S23) 



X^-^PP(^So)7rO,o('S23) 



ah 



:i7) 



Here h = tt^tt^ stands for the tt^tt^it^ production, and h = rjrj for -K^rirj. The /^-matrix terms 
which describe the prompt /o-production in the pp annihilation have the following form: 



-^p|(iSo)7rO,a("523) 



A 



(a) 

ppQ-So)-K° 



a) 



Ml - S23 



+ 0pp(i5o)7rO,a[OO] 



1 GeV^ + So \ / S23 - sa 



S23 + So 



^ S23 + Sao ' 



(18) 
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The parameters App(i5o)^o[00] and 0pp(i5o)7ro[OO] (or App(i5g)^[00] and 0pp(i5„)^[OO]) can be complex- 
valued, with different phases due to three-particle interactions. The matter is that in the final 
state interaction term (0) we take into account the leading (pole) singularities only. The next- 
to-leading singularities are accounted for effectively, by considering the vertices pp — >■ mesons 
as complex factors, for more detail see . 

Here the formulae are written for the production of two mesons in the (/ J^*" = 00^^)-state: 
this very wave has been fitted in the present analysis. In the final state of reactions under 
consideration, vector and tensor mesons are produced. As was said above, the parameters for 
vector and tensor resonances have not been fitted in this analysis; we have used those found 
]TT[] . The formulae for vector and tensor resonances used in fitting procedure are given in 



m 



Appendices A and B together with formulae for the annihilation from higher states. 



2.3 Quark-combinatorics rules for the decay couplings 

The decay couplings of the gg-meson and glueball to a pair of mesons are determined by the 
planar diagrams with gg-pairs produced by gluons: these diagrams provide the leading terms 
in the 1/N expansion 1^^, while non-planar diagrams give the next-to- leading contribution. 



The production of soft gg pairs by gluons violates flavour symmetry, with the following 
ratios of the production probabilities: 

uu : dd : ss = 1 : 1 : X . (19) 

Suppression parameter A for the production of strange quarks varies in the interval < A < 1. 
An estimate performed for the high energy collisions gives A = 0.4 — 0.6 in hadron decays 
it was evaluated as A ~ 0.8 [ESl. 



We impose for the decay couplings of bare states, the quark-combinatorics relations. 
The rules of quark combinatorics were previously suggested for the high energy hadron produc- 
tion 1 39 1 and then extended for hadronic J/ \l/ decays |^ . The quark combinatoric relations were 



used for the decay couplings of the scalar-isoscalar states in the analysis of the quark-gluonium 
content of resonances in and later on in a set of papers p, ^ |10|, [ll|, ^ We 



present the decay couplings for scalar mesons, which are a subject of our analysis, in Appendix 
C. 

The wave function of the /o-state is supposed to be a mixture of the quark-antiquark and 
gluonium componets as follows: 

gg cos a + gg sin a , (20) 

where the gg-state is a mixture of non-strange and strange quarks, nn = {uu + dd) / \/2 and ss: 

qq = nn cos (f + ss simp . (21) 

Using formulae given in Appendix C for the vertices gg —>■ vrvr, KK, rjr], rji]' together with 
analogous couplings for the transition gluonium— ^two-meson state, we obtain the following 
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coupling constants squared for the decays /o nn, KK, rji], rji]': 



2 ^ ( 9 , G 



2 



= 2 I f (sin V? + W ^ cos v?) + 



2' V2 V2 + A 

2 ^ I / COS^ /-- . -S/^x G / 2 \ ■ 2 

9m ~ n \ 9 \ — 7^ '^'^^ V9 + V A sm (y9 sm k3 ) -\ — (cos fc) + A sm c) 



,2 



1 rz- . ^ l-\ V 



g~^i = sin^ 9 cos^ 6 ( g {—^ cos (p — ^/X sin ip) + G 



V2 V^TxJ 

Here g = go cos a and G = Gosina, where go is universal constant for all nonet members and 
Go is universal decay constant for the gluonium state. The value gl^ is determined as a sum 
of couplings squared for the transitions to vr+vr" and vr'^Tr", when the identity factor for 7r°7r° is 
taken into account. Likewise, g'^j^ is the sum of coupling constants squared for the transitions 
to KK and K^K^. The angle stands for the mixing of nfi and ss components in the r] and 
Tj' mesons (see Appendix C). 

The quark combinatorics make it possible to perform the nonet classification of bare states. 
In doing that we refer to fo'^^'^^^s as pure states, either qq or glueball. For the qq states that 
means: 

(1) The angle difference between isoscalar nonet partners should be 90°: 

cp[ft^'-'\l)] - v[fr'\2)] = 90° ± 5° . (23) 

(2) Coupling constants go presented in Appendix C should be roughly equal to each other for 

all nonet partners: 

9o[f!^'^'"'\l)] ^ 9o[fi'^'"'\2)] ^ <?o[ar^^] ^ 9o[Kt'''^] ■ (24) 



(3) Decay couplings for bare gluonium should obey relations for glueball (see Appendix C). 



Conventional quark model requires exact coincidence of the couplings but the energy depen- 
dence of the decay loop diagram, B{s), may violate the coupling- const ant balance because of the 
mass splitting inside a nonet. The K-matrix coupling constant contains additional s-dependent 
factor as compared to the coupling of the N/D-amplitude |^: g'^{K) = g'^{N/D)/[l + B'{s)]. 
The factor [1 + B'{s)]~^ mostly affects the low-s region due to the threshold and left-hand 
side singularities of the partial amplitude. Therefore, the coupling constant equality is mostly 
violated for the lightest OO"*""*" nonet, l^Pg QQ- We allow for the members of this nonet 1 < 
g[fo{l)]/ g[fo{2)] < 1.3. For the 2^Po qq nonet members, we set the two-meson couplings to be 
equal both for isoscalar and isovector mesons. 
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3 Description of data and the results for the 00++-wave 

Table [l| demonstrates the data sets which have been fitted in the present analysis; here the 
number of fitted points for each reaction is also shown. 

For the description of the 00"*^+ wave in the mass region below 1900 MeV, five i^'-matrix 
poles are needed (a four-pole amplitude fails to describe well the data set under consideration). 
Accordingly, five bare states are introduced. As in the previous analysis we have found three 



solutions. All Solutions, which we denoted, following ||Tl|, Solutions I, II-l and II-2, are similar 



in the pricipal points to that found in [11 

The values for each reaction in every solution are given in Table 1. 

Figures 1-11 demonstrate the two- meson spectra and angle distributions for the reactions 
pp — > three mesons and np — >■ three mesons. The scattering amplitudes vrvr tttt, vrvr rjr], 



and so on, are restored on the same level of accuracy as in [|T^, and we do not show them. 



In Figs. 1, 2, 3 one can see the M^o^o, M^^o, M^^ distributions in the leactions pp{liquid H2) - 
^o^o^o ^0^^^ ^o^o^_ rjj^g Dalitz-plots for these reactions were the corner-stones for the K- 
matrix fits to the OO"*""*" wave at the early stage of our study P, ||: Figs 1, 2, 3 demonstrate the 
level of accuracy of the present fit. Figure 4 shows similar mass distributions in the reaction 
pp{liquid H2) 7r+7r°7r^. 

Figures 5 and 6 demonstrate the distributions for the reactions np{liquid D2) vr^vr^Tr", 
7r~7r~7r"'". The angle distributions stand for bands with the production of /2(1285) + /o(1300) 
(Fig. 5c) and p(1450) (Figs. 5d, 6d). In connection with the discussion of a possible resonance 
structure in the 1 = 2 channel, we show in Fig. 6c the angle distribution along the band 
M^-^- ~ 1300 MeV (recall that resonance states with 1 = 2 are not included into our fit). 

In Figs. 7 and 8, one can see the mass and angle distributions in the reactions pp{liquid H2) — 
KsKs'it'^, K~^K^tt^ (parameters for a'^"'^^^ and Kq'"^^^^ have been fixed, correspondingly, in ||lT| 
and 11 (JsT-matrix re-analysys of the SLAC-NAGO-CINC-INUS data || for Kn Kir). 



Angle distributions in Figs. 7d and 8d are selected for the demonstration, in connection with 
the discussed possibility for the ag-resonance to exist in the region ~ 1300 MeV (in the fit 
performed in |]lT|, the ao-state had not been found in the mass region ~ 1300 MeV). 



In Fig. 9, the distributions in the reactions pp{liquid H2) K lK^ tt~ {K lK~ tt^) , KsK^n^ 
are shown. The distributions in the channel KsK~7t~^ have not been fitted because of the 
acceptance-correction problem: the curve is the result of the fit carried out without taking 
account of this channel data; nevertheless, as is seen from the figure, the description is rather 
good. 

In Figs. 10 and 11, we show distributions in the reactions np{liquid D2) — > KsKs7r~ , KsK~tt^. 
The angle distributions (Figs. lOd, lid) are presented for the band where the production of 
Co (1300) can be seen; yet, we do not see any signal from ao(1300). 
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Table 1: 

List of reactions and values for the X-matrix solutions. 





Solution I 


Solution II- 1 


Solution II-2 


Number of 
points 




The Crystal 


Barrel data 


from liquid H2: 










pp — > 7r°7r°7r° 


1.30 


1.39 


1.32 


7110 


pp — > 7r°7777 


1.32 


1.30 


1.32 


3595 


pp — > TT^TT^T) 


1.22 


1.31 


1.34 


3475 


from gaseous H2: 










pp — > 7r°7r°7r° 


1.38 


1.35 


1.39 


4891 


pp —>■ TT^rjr] 


1.25 


1.24 


1.26 


1182 


pp — > 7r°7r°77 


1.17 


1.18 


1.20 


3631 


with charge pions, 










from liquid H2: 










pp TT+TT^TT" 


1.38 


1.42 


1.46 


1334 


from liquid D2: 










pn — > 7r°7r°7r~ 


1.46 


1.43 


1.48 


825 


pn — > 7r^7r"7r+ 


1.51 


1.57 


1.58 


823 


with kaons, 










from liquid H2'- 










pp KsKsTT^ 


1.08 


1.08 


1.10 


394 


pp K+K-n^ 


0.94 


0.91 


0.99 


521 


pp — > KlK^t:^ 


0.79 


0.81 


0.80 


737 


from liquid D2: 










pn KsKs7i~ 


1.73 


1.61 


1.69 


396 


pn — > KsK~TT^ 


1.33 


1.40 


1.22 


378 


CERN-Miinich data 










7r+7r~ TT'^TT" (all waves) 


1.82 


1.70 


1.68 


705 


GAMS data 










TTTT 7r°7r° (S-wave) 


1.15 


1.13 


1.28 


16 


TTTT — 777/ (S-wave) 


1.16 


1.35 


0.92 


16 


TTTT — > 7]?]' (S-wave) 


0.54 


1.76 


0.67 


8 


BNL data 










TTTT ^ A' A' (S-wave) 


1.15 


0.75 


1.10 


35 


SLAC-NAGO-CINC-INUS 










data 










K~7T^ — >• Ar~7r+ (S-wave) 


2.80 


2.80 


2.80 


46 
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3.1 Bare /g-states and resonances 

In the i^'-matrix analysis of the 00++-wave, five bare states have been found, see Tables 2, 3 
and 4. The bare states can be classified as nonet partners of the qq multiplets I^Pq and 2^Pq 
or a scalar glueball. The /^-matrix solutions give us two variants for the glueball definition: 
either it is a bare state with the mass near 1250 MeV, or it is located near 1600 MeV. 

After having imposed the constrains and (P^, we found the following variants for the 
nonet classification. 

Solution I: 

/o^^''"(700±100) and /o^'^'"^(1245±40) are I^Pq nonet partners with (^[/(^'^'■^(700)] = -70°±10° 
and V7[/(5'^'^<^(1245)] = 20° ± 10°. 

For members of the 2^Po nonet, there are two variants: 

1) either /(^'^'■^(1220 ± 40) and /(^'^''^(ITSO ± 30) are 2^Pq nonet partners, with (^[/($"^''^(1220)] = 
33°±8° and(^[/o^^'''=(1750)] = 60°±10°, while /o^'^'''^(1630±40) is the glueball, with ^[f}^^'%lQ^Qi)] = 
27° ± 10°; or 

2) /o^'^'"^(1630 ± 40) and /o^^''"(1750 ± 30) are 2^Po nonet partners, and /o^^'^^(1220 ± 40) is the 
glueball. 

Solution II-l: 

/o^^"''^ (670 ±100) and /o^^'''^(1215±40) are I^Pq nonet partners with v?[/o^^'^<=(670)] = -65°±10° 
and ip[f^^'%1215)] = 15° ± 10°; 

/o^''''^(1560 ± 40) and /o^^''*'(1830 ± 40) are 23Po nonet partners with (p[f^^'%15Q0)] = 15° ± 10° 

and iflf^'^'^lSSO)] = -80° ± 10°, 

/o^^''^(1220) is the glueball, ^[/o^^''<=(1220)] = 40° ± 10°. 

Solution II-2: 

/o^''''^(700±100) and /o^'^^'=(1220±40) are I^Pq nonet partners with <^[/(5'^'^(700)] = -70°±10° 
and ip[fl^^'%1220)] = 15° ± 10°. 

In this Solution there are two variants for the 2^Po nonet: 

1) either f^'^'^USO ± 30) and /o^^'^'=(1830 ± 40) are 23Po nonet partners with (^[/o^^'^<=(1230)] = 
43° ± 8° and ip[f^'''^{1830)] = -60° ± 10°, 

/o^^'"^(1560 ± 30) is the glueball, with iplf^'^'^lbQO)] = 15° ± 10°, or 

2) /(^^'■^(1560 ± 30) and /o^''''"(1830 ± 40) are nonet partners and /o^'''""(1230) is the glueball. 
Tables |^, | and H present parameters which correspond to these three Solutions. 

3.2 GAMS and E852 data on the tttt S'-wave production ir^p 
{7T7t)s n at moderate momenta transferred to nucleon, < |t| < 1.5 
(GeV/c)2 

The formulae of Section 2.2 allow one to describe simultaneously the tttt spectra in the reaction 
TT-p {7nT)sn at piab = 18 GeV/c pl and piab = 38 GeV/c [0. 
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Table 2: 

Masses, coupling constants (in GeV) and mixing angles (in degrees) for the /o^^-resonances 
for Solution I. The errors reflect the boundaries for a satisfactory description of the data. The 
sheet II is under the tttt and An cuts; the sheet IV is under the tttt, An, KK and 7777 cuts; the 
sheet V is under the nn. An, KK, rjrj and rjrj' cuts. 





Solution I 




a — 1 


a = 2 


a — 3 


a = 4 


a — 5 


M 


0.650l:Si° 


1.2451:°^° 


1.220l:°i° 


1.630l;°i° 


1.750l;°l° 


^(") 


0.940l:?go 




0.680l:g^° 


0.680l;°i° 


0.790l;°i° 


(a) 
^5 








0.9601;}°° 


0.9001;°^° 


0.2801;}°° 


ipaideg) 


-(72lfo) 


18.011 


331^ 


97+10 
^ ' -10 


-591}° 




(I = nn 


a = KK 


(I = ijii 


(I = ipi' 


(1 = In 


ha 


-0.0501;}°° 


0.250ti°° 

fba = 


0.4401;}°° 
6 = 2,3,4,5 


0.3201;}°° 


-0.5401;}°° 


sheet II 






Pole position 






sheet IV 




1.3061-- 
-z(0.147l:°^i) 


1.489l;°°« 
-^(0.05ll;°°i) 


1.4801;}- 
-.(1.030l;°?g) 




sheet V 


-i(0.072l:°}i) 
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Table 3: 

Masses, coupling constants (in GeV) and mixing angles (in degrees) for the /o^^-resonances for 
Solution II-l. The errors reflect the boundaries for a satisfactory description of the data. The 

sheet II is under the tttt and An cuts; the sheet IV is under the tttt, An, KK and 7777 cuts; the 
sheet V is under the nn. An, KK, rjrj and rjrj' cuts. 





Solution II- 1 






a — l 




a — 3 


q; = 4 


a — 5 


M 






i.ZiO_ 040 


1.2201:°^^ 


1.5601:°^° 


1.830l:°i° 


^(") 




0.9901:?^° 




0.6701:1^° 


0.5001:°^° 


0.4101:°^° 


(a) 
^5 










0.8701:1°° 


0.6001:1°° 


-0.850l:°i° 


ipaideg) 




-(66l?o) 




401}^ 


151°^ 


-8011° 






(I = nn 


a = KK 


a = ijii 


(I = ipi' 


(1 = In 


ha 






0.1001;°^° 

fba = 


0.3601:1°° 
6 = 2,3,4,5 


0.3201:1°° 


-0.3501:°°° 


sheet II 








Pole position 






sheet IV 






1.3201-- 
-z(0.130l:°^i) 


1.4851:°- 
-z(0.055l:°°i) 


1.5301:-° 

-^(o.9ool:Z) 




sheet V 


1.7851:1;- 

-^(0.1351:°?^) 
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Table 4: 

Masses, coupling constants (in GeV) and mixing angles (in degrees) for the /o^^-resonances for 
Solution II-2. The errors reflect the boundaries for a satisfactory description of the data. The 
sheet II is under the tttt and in cuts; the sheet IV is under the tttt, Att, KK and 7777 cuts; the 
sheet V is under the tttt, 47r, KK, rjr] and rjr]' cuts. 





Solution II-2 






a — 1 


a = 2 


a — 3 


a = 4 


a — 5 


M 




0.650l:Ji° 


1.220l:°^g 


1.230l:°i° 


1.560l:°i° 


1.8301:°!° 










0.4701:°^° 


0.4201:°^° 


0.420l:°t° 


(a) 
^5 










0.870li°° 


0.5601:°?° 


-0.7801:°?° 


(Pa{deg) 




-(68l?5) 


1411 


4311 


15112 


-00_io 






(I = TiTi 


a = KK 


a = ijii 


(I = ipl' 


a = Itv 


ha 






0.1001:1°° 

ha = 


0.2601:1°° 
6 = 2,3,4,5 


0.2601:1°° 


-o.i4ol:°i° 


sheet II 








Pole position 






sheet IV 






i.OZd_ 030 

-z(o.i7ol:^) 


1.4901:-;] 
-^(o.o6ol:°°i) 


1.4501:- 

-^(0.8001:1°°) 




sheet V 


1.7401-- 
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The leading vr and Oi trajectories have rather close intercepts and slopes [33|, so corre- 
sponding reggeon exchanges do not cause a change of t-distributions with energy. However, 
the experimental data definitely point to a break of similarity of the spectra: in Fig. 12 the 
difference of the M7r,r-distributions is shown for |t| ~ (0.3 — 0.4) (GeV/c)^ where one can trace 
differences in spectra related to E852 and GAMS experiments. This fact definetely points to a 
significant contribution of daughter trajectories. 

We have fitted to data and |]T^ under a variety of assumptions on the |t|-exchange 
structure; these versions of fitting are discussed in a separate publication |^. Here, in Figs. 13 
and 14, we demonstrate typical description of the data. Figure 13 represents the M^^^ spectra 
of E852 Collaboration for different |t|-intervals, and figure 14 shows the same for GAMS group 
(calculated spectra refer to Solution II-2). 

In Fig. 15, we show form factors for vr and ai exchanges, which at large |t| describe effectively 
the multi-reggeon interactions. We show the form factors which correspond to the so-called 
Orear behaviour of the scattering amplitude at large |t| (see also and references therein): 



97, 



exp 



+ A(t 



exp(-4"Vl^-"^^ 



(25) 



Another types of the behaviour of G^"-' (t) were analysed in [R^ . 



4 /o-resonances: masses, decay couplings and 
partial widths 

The resonance masses and decay couplings are not determined directly in the fitting procedure. 
To calculate these quantities one needs to perform analytical continuation of the ii'-matrix 
amplitude into lower half-plane of the complex plane s. One is allowed to do it, for the K- 
matrix amplitude correctly takes into account the threshold singularities related to the tttt, 
TTTTTTTr, KK, Tji], 1]!]' channels which are important in the 00"'""^-wave. 



4.1 Masses of resonances 

The complex masses of the resonances /o(980), /o(1300), /o(1500), /o(1200 — 1600) obtained 
in Solutions I, II-l and II-2 do not differ strongly. Solutions I and II are essentially different in 
the characteristics of the /o(1750) 

For the pole positions the following values, in MeV, have been found: 

Solution I : /o(980) ^ 1031 - i 32 (26) 



/o(980) 


1031 - 


i 32 


/o(1300) 


1306 - 


i 147 


/o(1500) 


1489 - 


i 51 


/o(1750) 


1732 - 


i 72 


/o(1200 - 1600) - 


1480 - 


i 1030 
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Solution II - 1 : /o(980) ^ 1020 - i 34 (27) 



/o(980) 


1020 


-i 34 


/o(1300) 


1320 


-i 130 


/o(1500) 


1485 


-i 55 


/o(lT50) 


1785 


-i 135 


/o(1200 - 1600) - 


1530 


930 



Solution II - 2 : /o(980) ^ 1020 - i 35 (2^ 



/o(980) 


1020 


-i 35 


/o(1300) 


1325 


-i 170 


/o(1500) 


1490 


60 


/o(1750) 


1740 


-z 160 


/o(1200 - 1600) - 


1450 


-z 800 



We see that Solutions I and II give different magnitudes for total width of the /o(1750). 



4.2 Decay couplings and partial widths 

We determine the couphng constants and partial decay widths using the following procedure. 
The 00"'""'"-amplitude for the transition a ^ b, 

Aa^b{s) , a,b = TT'K,KK,r]r],r]r]',7nnnr, (29) 

is considered as a function of the invariant energy squared s in the complex-s plane near the 
pole related to the resonance n. In the vicinity of the pole the amplitude reads: 

(n) (") 

Aa-,bis) = e'(^«"'+^^"') + non - pole terms . (30) 

Here /i„ is the resonance complex mass /i„ = M„ — ir„/2 ; g^^^ and f?^""* are the couplings for 
the transitions fo^a and /o b. In (pO]), the non-pole background terms are omitted. 

The decay coupling constants squared, g^'^ ^, are shown in Table 5 for a = vrvr, KK, rjr], rjrj', 
TTTTTrTT. The couplings are determined with the normalization of the amplitude used in Section 
2: for example, we write the vrvr scattering amplitude (^) as 

A^^^^^{s) = (r^oexp (2^5°) - l)/2ip^^{s) , 

where rj^ and 6q are the inelasticity parameter and phase shift for the 00++ Trvr-wave, reprec- 
tively. The coupling constants g^^ are found by calculating the residues of the amplitudes 
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Table 5: Coupling constants squared (in GeV^) of scalar-isoscalar resonances to hadronic chan- 
nels TTTT, KK, rjrj, rjr]' and ninnT for different K-matrix solutions. 



Pole position 


TTTT KK rji] rjT]' -k-k-r-k 


Solution 


/o(980) 
1031 - i32 
1020 - i35 
1020 - 235 


0.056 0.130 0.067 - 0.004 
0.051 0.115 0.051 - 0.003 
0.054 0.117 0.139 - 0.004 


I 

II- 1 
II-2 


/o(1300) 
1306 -il47 
1320 -il30 
1325 -il70 


0.036 0.009 0.006 0.004 0.093 
0.038 0.003 0.005 0.012 0.109 
0.053 0.003 0.007 0.013 0.226 


I 

II- 1 
II-2 


/o(1500) 

1489 - z51 
1485 - z55 

1490 - z60 


0.014 0.006 0.003 0.001 0.038 
0.020 0.007 0.004 0.003 0.049 
0.018 0.007 0.003 0.003 0.076 


I 

II- 1 
II-2 


/o(1750) 
1732 - z72 
1785 - il35 
1740 - a60 


0.013 0.062 0.002 0.032 0.002 
0.066 0.003 0.007 0.033 0.080 
0.089 0.002 0.009 0.035 0.168 


I 

II- 1 
II-2 


/o(1200 - 1600) 
1480 - ilOOO 
1530 - i900 
1450 - iSOO 


0.364 0.265 0.150 0.052 0.524 
0.325 0.235 0.086 0.015 0.474 
0.179 0.204 0.046 0.005 0.686 


I 

II- 1 
II-2 



TTTT — >■ TTTT, KK, rjT], TjT]' , TTTTTrTT. Also we check the factorization property for the pole terms by 
calculating residues for other reactions, such as KK —>■ KK, and so on. 

The resonance /o(980) is located near the strong KK threshold, therefore two poles are 
related to /o(980): for example, in Solution II-2 the nearest one is on the 3rd sheet (under vrvr 
and TTTTTTTr cuts, at M ~ 1020 — i35 MeV) and a remote pole on the 4th sheet (under vrvr, vrTTTTTr 
and KK cuts, at M ~ 940 — ^240 MeV). Coupling constants for /o(980) are determined as 
residues of the nearest pole which is located on the 3rd sheet. 

The partial width for the decay /o ^ a is determined as a product of the coupling constant 
squared, g^^ ^, and phase space, Pa(s), averaged over the resonance density: 

The resonance density factor, [{Re /i^ — s)^ + (Im/i^)^] ^, guarantees rapid convergence of the 
integral (0). The normalization constant C„ is determined by the requirement that the sum 
of all hadronic partial widths is equal to the total width of the resonance: 

r(n) = Er,(n). (32) 
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Table 6: Partial widths of scalar-isoscalar resonances (in MeV) in hadronic channels tttt, KK, 



rjrj, rjrj' and irinTTC for different i^-matrix solutions. 





TTTT 


KK 


r]r] 


rjrj' 


TTTT TTTT 


Ttot/2 


Solution 


/o(980) 


52 


10 


— 


— 


2 


32 


I 




57 


9 


- 


- 


2 


34 


II-l 




58 


10 


- 


- 


2 


35 


II-2 


/o(1300) 


77 


14 


7 


— 


196 


147 


I 




66 


4 


4 


- 


186 


130 


II-l 




63 


3 


5 




269 


170 


II-2 


/o(1500) 


23 


6 


3 


0.0 


70 


51 


I 




25 


7 


3 


0.1 


75 


55 


II-l 




22 


6 


2 


0.1 


88 


59 


II-2 


/o(ir50) 


31 


96 


5 


7 


5 


72 


I 




114 


3 


8 


7 


138 


135 


II-l 




103 


1 


8 


4 


204 


160 


II-2 



In Table 6 we show the values of partial widths for the resonances /o(980), /o(1300), /o(1500), 
/o(1750). Partial widths for /o(1300), /o(1500), /o(1750) are calculated within standard formu- 
lae for the Breit-Wigner resonances (pTI). 



Because of a strong KK threshold near the pole, the resonance term for /o(980) should 
be used not as the Breit-Wigner amplitude but in a more complicated form. Instead of the 
Breit-Wigner pole term i?^"*-* = g'^^'g'if' / ilJ-n — s) entering equation (0), the following resonance 
amplitudes can be written for the tttt —>■ tttt and KK —>■ KK transitions near /o(980) |^ : 



R 



(tttt, tttt) 
/o(980) 



s — 4m|- 



1 

'd 



(33) 



where 



R 



{KK,KK) 
/o{980) 



Gl^ + z{2GGj,Kf + f 



.rrin 



s)) 



1 

D 



D 



niQ ~ s — iG — i 



s — 4m|- 



Glji + z{2GGj,jif + f{ml-s)) 



(34) 



Here mo is the input mass of /o(980); G and G^^ are coupling constants to pion and kaon 
channels. The dimensionless constant / stands for the prompt transition KK tttt: the value 
//mo is the "transition length" which is analogous to the scattering length of the low-energy 
hadronic interaction. The constants mo, G, G^Ri f ^^e parameters which are to be chosen 
to reproduce the /o(980) charactristics. Equations (^) and (|34D are written for s > 4:mj^, at 
s < 4m|- one needs to replace 



4m|- 



4m|^ 



s. 
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In two sets of parameters were obtained, with sufficiently correct values for the /o(980) 
pole position and couplings. They are equal (in GeV) to: 



Solution A: mo = 1.000 , / = 0.516 , G = 0.386 , Gkr = 0.447 , 
Solution B : mo = 0.952 , / = -0.478 , G = 0.257 , G^r = 0-388 . 



(35) 



Partial widths of /o(980) are calculated with the expression similar to (BH 
placement of the integrand denominator as follows: 



(Re/i^ — s)^ + (Imyu: 



2\2 



\D\ 



with the re- 



(36) 



For both sets of parameters ( pSj) the calculated partial widths are close to each other. The 
values of partial widths for /o(980) averaged over Solutions A and B are presented in Table 6. 

The total hadron width of /o(980) is defined in the same way as for the other /o-mesons, 
namely, by using the position of pole in the complex-M plane: the imaginary part of the mass 
is equal to a half-width of the resonance. For the Breit-Wigner resonance this definition is in 
accordance with what is observed from resonance spectrum (provided there is no interference 
with the background). If the resonance is located in the vicinity of a strong threshold, the 
visible resonance width differs significantly from that given by pole position. 

The results of our calculations of partial decay widths are presented below, the magnitudes 
are in MeV. In the brackets we demonstrate the values obtained in on the basis of previous 
analysis ||Tl[| : 



Resonance 
/o(980) : 

/o(1300) : 

/o(1500) : 

/o(1750) : 



r 

55 ±5 
(64 ±8) 
66± 10 
(46 ± 12) 

23 ±5 
(37 ±2) 

30 ±5 
105 ± 10 

(741^^) 



^KK 
10 ± 1 

(12 ±1) 
6±4 
(5 ±3) 
6±2 
(7 ±2) 
100 ± 10 
2± 1 



vv 



5±2 
(4 ±2) 
3± 1 

(4±1) 
5±3 
8± 1 

(7±1) 



0.1 ±0.1 
(0.2 ±0.1) 
7±3 
5±3 
(3±1) 



2± 1 
(4 ±2) 
230 ± 50 
(171 ±10) 
80± 10 
(79 ± 6) 

5±5 
170 ±40 
(91: 



+30 N 
-60 J 



Ttot/2 

34 ±5 
(40 ± 5) 
150 ± 20 
(113 ±10) 

55 ±5 
(62 ± 3) 
70 ±20 
150 ± 25 

(931^°) 



Solution 
IJI 



11, 49 



IJI 



11, 49 



I,II 



IT], Si 
I 

II 



11, 



(37) 



In Solutions I and II, partial widths for /o(980), /o(1300), /o(1500) are agreeably coincide but 
for /o(1750) the decay characteristics are different. So we present separately the /o(1750) partial 
widths for both Solutions, I and II. 



Classification of scalar resonances 



The X-matrix approach works ingenuously with the bare states thus directly providing the 
gg-classification for these states. Concerning the resonances, the classifications of the and 
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Kq mesons in terms of bare states and real resonances are similar, while /o mesons give us 
an opposite case: the matter is that there is an exotic state, gluonium, in the region under 
consideration, whose mixing with qq state is not forbidden by the 1/N expansion rules (for 
more detail see 0, 0). 

In the present Section, on the basis of the values for couplings to channels tttt, KK, tjt], rjr]' 
calculated in Section 4, we analyse the quark-gluonium content of resonances /o(980), /o(1300), 
/o(1500), /o(1750) and the broad state /o(1200 - 1600). 



5.1 Overlapping of the /o-resonances in the mass region 1200—1700 
MeV: accumulation of widths of the qq states by the glueball 

The existence of the broad resonance is not an eventual phenomenon: it originated due to the 
mixing of states in the decay processes, namely, transitions /o(^i) — ^ real mesons /o(^2)- 
These transitions result in a specific phenomenon, that is, when several resonances overlap, 
one of them accumulates the widths of neighbouring resonances and transforms into the broad 
state. 

This phenomenon had been observed in H, p| for scalar-isoscalar states, and the following 



scheme has been suggested in Q: the broad state /o(1200 — 1600) is the descendant of 
the pure glueball which being in the neighbourhood of qq states accumulated their widths and 
transformed into the mixture of gluonium and qq states. In |^2[, this idea had been modelled for 
four resonances /o(1300), /o(1500), /o(1200 — 1600) and /o(1750), by using the language of the 
quark-antiquark and two-gluon states, qq and gg: the decay processes were considered as the 
transitions /o — >■ qq,gg, correspondingly; the same processes are responsible for the mixing of 
resonances. In this model, the gluonium component is mainly dispersed over three resonances, 
/o(1300), /o(1500), /o(1200 — 1600), so every state is a mixture of qq and gg components, with 
roughly equal percentage of gluonium (about 30-40%). 

Accumulation of widths of overlapping resonances by one of them is a well-known effect 
in nuclear physics ||5l|, ^ In meson physics this phenomenon can play rather important 



role, in particular for exotic states which are beyond the qq systematics. Indeed, being among 
qq resonances, the exotic state creates a group of overlapping resonances. The exotic state, 
which is not orthogonal to its neighbours, after accumulating the "excess" of widths turns into 
the broad one. This broad resonance should be accompanied by narrow states which are the 
descendants of states from which the widths have been taken off. In this way, the existence of 
a broad resonance accompanied by narrow ones may be a signature of exotics. This possibility, 
in context of searching for exotic states, was discussed in [10, 



The broad state may be one of the components which form the confinement barrier: the 
broad states, after accumulating the widths of neighbouring resonances, play for these latter the 
role of locking states. Evaluation of the mean radii squared of the broad state /o(1200 — 1600) 
and its neighbouring resonances argues in favour of this idea, for the radius of /o(1200 — 1600) is 
significantly larger than that of /o(980) and /o(1300) [l^, |5^, thus making /o(1200 — 1600) to be 
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the locking state. The data on the t-dependence in the reaction n p {tth) sn at Pbeam = 18 
(GeV/c)^ ||19| also support this idea, for the signal from the broad state (background in the 



spectra) falls with the increase of \t\ much quicker than that of the /o(980) and /o(1300). 

The /^-matrix solutions obtained here (Section 3.1) give us two variants for the transfor- 
mation of bare states into real resonances coming after the onset of the decay channels. They 
are as follows: 



/o^»-(700)±100) ^ /o(980) , 

/o^'^^^(1220±40) ^ /o(1300) , 

/o^'^^^(1230 ± 40) ^ broad state /o(1200 - 1600) , (38) 

/o*'^'^^(1580±40) /o(1500) , 

/o^'^'^^(1800±40) ^ /o(1750) . 



and 



/o'«'-^(700)±100) ^ /o(980) , 

/o^"^^(1220±40) ^ /o(1300) , 

/o^"^^(1230±40) ^ /o(1500) , (39) 

/o^"'^'=(1580 ± 40) ^ broad state /o(1200 - 1600) , 

/o^'^'^^(1800±40) ^ /o(1750) . 

The evolution of bare states into real resonances is illustrated by Fig. 16: the shifts of amplitude 
poles in the complex-M plane correspond to a gradual onset of the decay channels. Technically 
it is done by replacing the phase spaces pa for a = tttt, TrvrTTTr, KK, rjr], rji]' in the i^'-matrix 
amplitude as follows: pa — »■ ^Pa, where the parameter ^ runs in the interval < ^ < 1. At 
^ ^ one has bare states, while the limit ^ ^ 1 gives us the positions of real resonances. 

5.2 Hadronic decays and evaluation of the quark-gluonium content 
of scalar-isoscalar resonances 

Here, on the basis of the quark combinatorics for the decay coupling constants, we analyse 
quark-gluonium content of resonances /o(980), /o(1300), /o(1500), /o(1750) and the broad state 
/o(1200 - 1600). 

The comparison of Tables 2, 3, 4 with Table 5 demonstrates a strong change of couplings 
during the evolution from bare states to real resonances; the same type of strong deviation of 



couplings has been observed before, for previous i^-matrix solutions, see ||5^. Note that the 
change occurs not only in absolute values of couplings but in relative magnitudes as well that 
means the change of the quark-antiquark content of states in the evolution caused by the decay 
onset. 

Let us look at the proportion of ss, nn and gluonium components given by quark-combinatorics 



relations (p2|) for the studied resonances. The coupling constants squared for /o mr, KK^ rjr], 
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rjr]' in (^2]) are expressed as a sum of two terms which correspond to transitions of quarkonium 
and gluonium components into two pseudoscalar mesons. 

First, let us find the mean value of mixing angle, (if), for the nn/ss components in the 
intermediate state. Because of the two-stage mechanism of the gluonium decay gluonium — >■ 
quark — antiquark pair — > two mesons (see ^ and Appendix C), we determine (ip) in the 
intermediate state as follows: 

/o gluonium + qq nn cos{ip) + sssm{ip) two mesons. (40) 

So, {{p) is the angle corresponding to the coupling constants squared (|22|) at G = 0. 

Based on the values of couplings given in Table 5, we have found the following values of (ip) 
and A. 
Solution I: 



/o(980) 
/o(1300) 

Broad state /o(1200 - 1600) 
/o(1500) 
/o(1750) 



(V?) ~ -68° , A ^0.5 -1.0 

((^) ~ (-3°) - 4° , A ^0.5 -0.9 

{if) ~ 27° , A ^ 0.54 , 

(^)~12°-19°, A ~ 0.5 -1.0 

{ip) ~ -72° , A ~ 0.5 - 0.7 



(41) 



Solution II- 1: 



/o(980) 
/o(1300) 

Broad state /o(1200 - 1600) 
/o(1500) 
/o(ir50) 



{(p) ~ -67° , A ~ 0.5 - 1.0 , 

((^) ~ (-21°) - (-10°) , A ~ 0.5 -1.0, 

{ip) ^ 28° , A ~ 0.55 , 

{ip) ~ 0° - 11° , A ~ 0.5 - 1.0 , 

{ip) ~ -37° , A ~ 0.45 - 0.55 , 



(42) 



Solution II-2: 



/o(980) 
/o(1300) 

Broad state /o(1200 - 1600) 
/o(1500) 
/o(1750) 



((^) ~ -67° , A ^ 0.6 - 1.0 

(V^) ~ (-16°) - (-13°) , A ~ 0.5 -0.6 

{ip) ^ 33° , A ~ 0.85 , 

((^)~2°-ir, A ^0.6 -1.0 

{ip) ~ -18° , A ^ 0.5 , 



(43) 



Note that the relations (^21) provide one more solution for /o(980), with {ip) ~ 40°. However 
this solution contradicts to v9[/q"'^'^(700)] ^ —67° obtained within the K-mathx fit (Tables 3, 4, 
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5): for details see [^], where the transformation of bare poles into resonances after the onset 
of the decay channels, was traced. 

When the resonance /o is considered as a system gluonium + gg, where qq = uu cos ip + 
sssmip, the coupling constants given in Table 5 determine as a function of the ratio G/ g. The 
results of fitting to coupling constants for comparatively narrow resonances /o(980), /o(1300), 
/o(1500), /o(1750) and the broad state /o(1200 — 1600) are presented in Fig. 17. In Fig. 17a,c,e, 
the curves present as a function of G/ with different A's from the intervals given by Eqs. 
(41)-(43) for the resonances /o(980), /o(1300), /o(1500), /o(1750): a bunch of curves provides 
us the values ((^9, G /g) which describe well the couplings squared of Table 5. One can see that 
the correlation curves for all resonances /o(980), /o(1300), /o(1500), /o(1750) are, with a good 
accuracy, the straight band. As is stressed below, such a behaviour of correlation curves is a 
signature of the qq origin of the resonances. 

The magnitudes g^ and G^ are proportional to Wqq and Wgiuonium which are the probabilities 
for quark and gluonium components to be present in the considered resonance: 
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_ /-i2 Ti^ 
^ ^ gluonium' ^ gluonium 



(44) 



According to the rules of 1/A^ expansion [^, the coupling constants, g'^^ and are of 



the same order (for more detail see [^, ^), therefore we accept as a rough estimation: 

G^/g^ ~ Wgiuomum/W,q . (45) 



Varying G/g in the interval —0.8 < G/g < 0.8 corresponds to a plausible admixture of the 
gluonium component up to 40%, Wgiuonium ^ 0.40. 

Figures 17a, c,e provide the following intervals ip in the resonances /o(980), /o(1300), /o(1500), 
/o(1750), after the admixture of the gluonium component. 

Solution I: 



WgluoniumiWm ^15% : 
WgluoniumMl^m ^ 30% 

W,;™™[/o(1500)] <30% 
Vr,,„™[/o(1750)] <30% 



-93° < </^[/o(980)] < -42°, 
-25° < <^[/o(1300)] < 25° , 
-2° < <^[/o(1500)] < 25° , 
-112°<<^[/o(1750)] <-32° 



(46) 



Solution II- 1: 



WgiuomumWm] ^ 15% 
Wgluonium[fo{l^m ^ 30% 
WgiuoniumMl^m ^ 30% 

W^5z™m[/o(1750)] < 30% 



-92° < <^[/o(980)] < -42°, 
-54° < </^[/o(1300)] < 13° , 
-19° < </p[/o(1500)] < 21° , 
-73°<<^[/o(1750)] <-2° . 



(47) 
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Solution II-2: 



WgiuomumlWm < 15% : 
WgluoniumlfoiliOO)] < 30% 

W^gz™™[/o(1500)] < 30% 

Wgluonium[fo{l750)] < 30% 



-90° < vlWSO)] < -43°, 
-42° < v^[/o(1300)] < 10° , 
-18° < ¥^[/o(1500)] < 23° , 
-46° < ^[/o(1750)] < 7° . 



The broad state /o(1200 — 1600) demonstrates another type of the {ip, G/g) correlation, for 
all Solutions I, II-l, II-2, when the constraints (^) are imposed on the coupling constants of 
Table 5 (see Fig. 17b, d,f). The correlation curves form a typical cross which, as is seen from 
the explanation given below, is the glueball signature. Indeed, the glueball can mix with the qq 
state because of the transition glueball {qq)giuebau — ^ glueball. The state {qq)giuebaii is close 
to the the SU(3) flavour singlet (qq) singlet = (uu + dd + ss)/a/3 but is slightly different. In the 
transition glueball {qq)giuebau glueball the production of the ss pair can be suppressed, 
and it looks reasonable to assume that this suppression is of the same order as it is in meson 
qq decays, where the the new qq pair is created by the gluon field: uu : dd : ss = 1 : 1 : X, with 
A ~ 0.5 — 0.8. Then {qq)giuebau = {uu + dd + \/X ss) /\/2 + X In terms of the nn and ss 
states, 

(qq) glueball = UH COS ^glueball + SS siu (figlueball , (49) 

where fgiuebaii = tg~-^y^A/2 ^ 27° — 33° for A ~ 0.5 — 0.8. The glueball descendant is a mixture 
of the gluonium (gg) and quarkonium {qq)giuebaii components: 

gg cos 7 + {qq)giuebaii sin 7 (50) 

Since the ratios of couplings for the transitions [gg) —>■ nn, KK, rjrj, rjrj' and {qq)giuebaii — *■ tttt, 
KK, rjr], rjr]' are quite the same, one cannot find out the mixing angle 7. But this very property 
- the similarity of ratios of coupling constants for the gluonium and quarkonium component - 
gives rise to special form of correlation curves on the {(p, g/G) -plot, the cross in Fig. 17b, d,f. 
Vertical part of the line means that the glueball descendant may have any noticeable proportion 
of the (qq)giuebaii state. Horizontal part of the curve corresponds to the dominantly gluonium 
component. 

The couplings from Table 5 for /o(1200 — 1600) allow certain deviation in the qq component 
from the state {qq) giuebaii- the couplings may be described by the values {g/G) and ip belonging 
to hyperboles of the type shown in Fig. 17b, d,f by dashed and dot-dashed lines. We have a 
bunch of such hyperboles which go off the centre of the cross, when A deviates from its central 
value (A = 0.55 for Solution I, A = 0.54 for Solution II-l and A = 0.85 for Solution II-2). 

An appearence of the gluonium cross on the {ip, g/G) -plot, when the formula (|22|) is used 
to determine coupling constants, is the glueball (or glueball descendant) signature in case of a 
strong mixing of the gluonium and quarkonium components. On the contrary, the absence of 
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the gluonium cross in the correlation curves should point to a qq origin of the considered /o 
state. 

Therefore, the broad state /o(1200 — 1600) in all Solutions is the glueball descendant, for 
the transition couplings /o(1200 — 1600) tttt, KK, rjri, rjrj' point directly to this fact. 

The states /o(980), /o(1300), /o(1500), /o(1750) cannot pretend to be the glueball descen- 
dants. 

The states /o(1300) and /o(1500) are dominantly the nn states. Though, in Solution II-l, 
the resonance /o(1300) is allowed to have a large ss component: with a noticeable admixture 
of the gluonium component, the angle v9[/o(1300)] may achieve the value ^ —50°. As to 
/o(1500), v9[/o(1300)] may reach 25° at G/g ~ —0.6 (Solution I). The description of couplings 
squared, g^, in this case owes a strong destructive interference of the decay amplitudes {qq) 
two pseudoscalars and (gg) — > two pseudoscalars. So in this case we should not be tempted by 
the proximity of (/)[/o(1500)] and (pgiuebaii to identify /o(1500) as the glueball or its descendant. 



5.3 Systematics of scalar states on the (n, M^)-plane 



The systematics of resonances carried out in demonstates that all resonances can be plotted 
on linear trajectories at the (n, M^)-plane, = Mq + — l)yU^, with a universal slope /x^ ~ 1.3 
GeV2. 

This empirical property of qq states may serve as an additional signature for the qq origin 
of resonances /o(980), /o(1300), /o(1500), /o(1750). These resonances fit well to linear trajec- 
tories, with the slope /x^ ~ 1.3 GeV^. The figure 18a demonstrates the (n, M^)-trajectories 
for resonance states with 00"*"+, 10++ and |0++, if the /o(1200 — 1600) is accepted to be of 
the glueball origin (recall that a doubling of the /o-trajectories occurs, due to the existence of 
two components, nn and ss). Similar trajectories for bare states are shown in Fig. 18b, if the 
/q"'''^(1580) is the gluonium. The trajectory slopes for real and bare states almost coinside. 

In Fig. 18, we present the variant where /q°^*^(1580) is the gluonium and the broad state 
/o(1200 — 1600) is its descendant; in this way it is natural that there is no room for these states 
to be on linear trajectories on the [n, M^)-plane. 

For bare states of Solution II-l, the linearity of trajectories is completely broken, see Section 
3.1, because the linearity exists for bare states only when /q"^*^(1580) is the glueball. This is 
an argument against Solution II-l as physical solution. 



6 Conclusion 

The i^'-matrix analysis of the 00++-wave based on the use of the spectra tttt, KK, rji], rjr]' in a 
broad variety of reactions (Table 1) provided us with three solutions: Solution I and Solutions II- 
1, II-2 (Tables 2,3,4). All these Solutions, despite a significant increase of the used experimental 



information, occurred to be similar to Solutions obtained in the previous analysis [11]. In all 
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Solutions, five poles have been found for the 00+"*" amplitude at complex masses in the studied 
mass region, 280 < M < 1900 MeV; they correspond to five scalar-isoscalar resonances, four 
of them being comparatively narrow resonances, /o(980), /o(1300), /o(1500), /o(1750), and a 
broad state /o(1200 - 1600). 

All Solution prove that comparatively narrow resonances /o(980), /o(1300), /o(1500), /o(1750) 
are of the qq origin. Both the ratios of the decay couplings for /q — > tttt, KK, rjr), rjr)', which 
result in linear behaviour of the correlation function {ip, G/g) (see Section 5.2), and the creation 
of linear trajectories in the (n, M^)-plane by resonances (Section 5.3) point to this fact. 

The broad state /o(1200 — 1600) has a gluonium origin: this is testified by 

(i) a specific behaviour of the correlation curve — glueball cross — in the {ip,g/G) plot (Section 

5.2) , and 

(ii) the absence of a room for the broad state on linear trajectories in the (n, M^)-plane (Section 

5.3) . 

Solutions II-l and II-2 provide us with proximate properties for all resonances /o(980), 
/o(1300), /o(1500), /o(1750), and a broad state /o(1200- 1600) but they differ by the classifica- 
tion of bare states, i.e. by the "formation history" of these resonances. In Solution II-l, a pure 
glueball state is located around 1200 — 1250 MeV, that results in the fact that bare gg-states 
do not form linear (n, M^)-trajectories. As to Solutions I and II-2, the pure gluonium state can 
be located around 1500 — 1600 MeV, so the bare qq states are rehably set on linear trajectories. 
The criterion of linearity for bare qq states allows us to consider as physical Solutions I and 
II-2. 

Solution I differs from Solution II by the content of /o(1750) which is dominantly ss state 
for Solution I and dominantly nn state for Solution II. 

The i^-matrix amplitudes reconstructed in our analysis do not contain the pole associated 
with a comparatively narrow sigma-meson. However, one should keep in mind that the K- 
matrix technique does not allow one to restore analytical amplitude in the region neighbouring 
left-hand cuts; this region is evaluated as Res < 4m^ or, in terms of the invariant mass M, 
as (ReM)^ — (ImAf)^ < 4m^. Because of that our X-matrix analysis is unable to provide 
definite conclusion about the existence of the a-meson with a large width. Nevertheless, let 
us stress that in some variants of the fit the restored i^T-matrix amplitudes have poles which 
might be considered as a light and broad cr-meson. For example, in Solution I such a pole is 
located at ~ 0.25 — i 1.00 GeV^, while we do not see similar pole around Im ~ 0.5 — 1.5 
GeV^ in Solution II. But, underline once again, this region cannot be treated as reliable in 
reconstructing analytical amplitude. 

We have analysed rich experimental information, yet it allowed us qualitative evaluation 
only of the percentage of the nn, ss and gluonium component in studied /o-mesons. In this 
sense, the information provided by hadronic decays of /o-mesons is exausted: the increase of 
data accuracy or of a number of studied reactions cannot help us to attain a considerable 
progress in understanding of /o-mesons in the region under investigation. A qualitative and 
precise reconstruction of the content of /o-mesons can be done after analysing a broad variety 
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of non-hadronic reactions. It looks like the study of reactions such as 77 — >■ tttt, KK, rjri, rjri' 
as well as J/\E' —>■ 'y'irn,'yKK,'jriri,'yrirj' (on the ground of relevant statistics) could provide 
the reliable and final magnitudes for the nn, ss and gluonium component in /o-mesons, for 
the production of hadrons in these reactions is initiated by the qq and gluonium components, 
correspondingly. 
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Appendix A: Amplitudes for the partial waves 02++, 10++, 
12++ 

Here we present the amplitudes for the waves which have not been fitted in the present analysis. 



The parameters of the resonances in these waves were found in [jTT[, and they are used as fixed 
valyes. 

Scattering amplitudes 

Scattering amplitudes for the partial waves 02++, 10++, 12++ are written as: 

i(") = i^(^-^)(i-^pK("))~' . (51) 



Isoscalar-tensor, 02++, partial wave 

The D-wave interaction in isoscalar sector is parametrized by the 4x4 i^-matrix where 1 = vrvr, 
2 = KK, 3 = 7]'!] and 4 = multimeson states: 

= D.(s) ( E + ^l.2) lGeV^ + -^A ^^(,) , (52) 

' ^ - S S + S2 



Factor -Da(s) stands for the D-wave centrifugal barrier. We take this factor in the following 
form: 

^«(^) = ^T^' « = 1,2,3, (53) 
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where ka = y s/A — m\ is the momentum of the decaying meson in the centre-of-mass frame of 
the resonance. For the multi-meson decay, the factor D/^^s) is taken to be 1. The used phase 
space factors are the same as those for the isoscalar S'-wave channel. 

Isovector-scalar, 10++, and isovector-tensor, 12++, partial waves 

For the amplitude in the isovector-scalar and isovector-tensor channels, we use the 4x4 K- 
matrix with 1 = vrr], 2 = KK, 3 = vri]' and 4 = multi-meson states: 

Al^W-^.wfef^ + /.i^^^]«.W. (54) 
' ^ Mf' - s s + si I 



Here J = 0,2; the factors -Da(s) are equal to 1 for the 10++ amplitude, while for the D-wave 
partial amplitude the factor -Da(s) is taken in the form: 

Da{s) = 1 . (55) 



Three-meson production amplitudes 

The partial waves 02++, 10++, 12++ are taken into account in the three-meson production 
processes. 

Invariant production amplitude A^^i j^2s+i£,^i b)(23) for the transition pp (^'^+^Lj) — > 1 + (2 + 
3)ij where the indices Ij refer the isospin and total angular momentum of the mesons 6 = 2 + 3 
{Ij = 02 10, 12) reads: 



^p^(2S+iLj)i,b(23) — X! -^pi (2S+iL_^)i,fe('523) 



ah 



(56) 



Here 

-^p|(2S+lLj)i,a('523) 



+ <^pp(2S+lLj)l,a[-^j] 



1 GeV^ + So 

S23 + So 



Da{s23), (57) 



and parameters ^^^(^2s+iij)i[Ij] 4>pp ('2s+iLj)i,a[Ij] Kiay be complex magnitudes with different 
phases due to three-particle interactions. 

Equations (|5^) and (|57D give us invariant parts of the final-state interaction block. The 
factors related to the angular-momentum expansion are presented in Appendix B. 
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Appendix B: Moment-operator expansion for the pro- 
cesses pp three mesons 

We analyse the processes pp — *■ three mesons by using the technique of moment-operator 
expansion; in this way our work is grounded upon the papers ^ |TU|, |TT[ . Below the necessary 
formulae are given for the reaction pp — > three mesons in the liquid and gas, when the pp 
annihilation is going from the lowest waves, S and P and the meson pair is in the S-, P- and 



D-states. In [^], this technique is presented in its general form, one may address this paper 
for more detail. 

The amplitude for the cascade transition pp resonance + meson three mesons has 
the following structure: 

The operator '4'{~(l2)Q^^^I^2^}^j -it^j'^^^^) refers to the pp state with the total angular momentum 
J: here S is the total spin of fermions, S* = 0, 1, and L is the angular momentum; qi and q2 
are the momenta of fermions and q = (qi — q'2)/2. The operator Q|j{'^2...^j_^^j(/i;i2, ^3) stands 
for the three-meson operator with the production of particles 1 and 2 in the resonance j-wave; 
this operator depends on relative momentum, ki2, of mesons 1 and 2, ku = (^1 — ^2)/2, 
and k^. Partial-wave amplitude is a function of the invariant energy squared s = P^, where 
P = ki + k2 + ks = qi + g2, and S12 = (fci + ^2)^- 

Angular-momentum operator for two particles, Xj^^jj^^ 

Let us introduce the operator of angular momentum of two particles, Xj^^jj^^^ ^^__^^^(q), which 
is constructed by using relative momentum of mesons in the space orthogonal to the total 
momentum P: 

P P 

(If, = (lu9ui,i 9vf, = 9vf, —■ (59) 

In the centre-of-mass system, where P = {Po,P) = (^5,0), the vector q^ is space-like: q^ = 
(0, g). We determine the operator X^^^j^^ ^^ _^^^{q) to be symmetrical and traceless. It is easy 
to construct it for the lowest values of L, L = 0, 1, 2: 

X''' = 1 , xl!^ = it , 41 = ? ( lUi - \ ilai^ , (60) 



Correspondingly, the generalization of Xj^^'' for L > 1 reads: 

X(^) = a^Z(^-i) (61) 
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2 ^ 



— Y X(^-i) 

It is seen that the operator Xj^^j^^^ ^^^^^^{q) constructed in accordance with (|61|) is symmetrical, 

X(^) = X^^^ (62) 

and it works in the space orthogonal to P: 

p ) „ ,^ = . (63) 
The angular-momentum operator X^^^ is traceless over any two indices: 

tlfJ-ifJ-j^'-lll...fli...llj...llL tlfJ.iflj^'-tll...fli...llj...llL \^^J 



Three-particle production amplitude with resonance in the interme- 
diate state 

The moment-operator for the production of three spinless particles in the cascade process 
resonance(j^'") + meson, can be written in terms of operators X^^^ To be definite, we 
consider the simplest case of vector resonance (j = 1) and give a generalization for j = 2. 

For the two-stage reaction of the type pp VP PPP, we denote the particle momenta 
of pseudoscalars (P) in the final states as ^1,^2,^3; the vector resonance (V) is produced in 
the channel 1 + 2 with the total momentum p = ki + k2 and relative momentum for the decay 
products ki2 = — k2). Then the final-state moment operator in (p8|) reads: 

Q£''.'L-i..(^i2, ks) = AW=J+')(s, s^2)Xlit'J,M)X^:\ki,) (65) 

+A^''^-'-'\s,s,2)Xj^^^_^iki)Xi^iki^, 
where ki2^ is orthogonal to p = ki + k2, 

kut^ = (a,.,!' - ) ^i2m' • (66) 
Here p^ = S12, and the momentum k^ is orthogonal to the total momentum P = ki + k2 + k^, 

ki = {9,,' - ^) ^3.' . (67) 
One has for the state with / = L: 

AW='\s, Su)Xj,il^^_^^{ki)xj^\ki,)P,e^p,,^ , (68) 
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where e^/s-y^j is the totally antisymmetrical four-tensor. 

The generalization for a higher resonance is obvious: for the tensor resonance, j = 2, in 
the reaction pp —>■ TS — * PPP, instead of two terms in (|65D one has three terms for the 
partial-wave amplitude: 

AW=J^')^s, s,,)X\f^,%^^^^{ki)X^SM2) (69) 

Moment operator for proton-antiproton system, ^2)Q^f'^/tjV(9i)) 

The partial-wave vertex for this system with the total angular momentum J, orbital momentum 
L and total spin S is determined by bilinear form Q'2)Q^^'/'2.'^./i7_iM/^(^i) where ('?/'(— ^2) and 
tpiQi)) are bispinors and Q'^^^l^^'^},^j_^^j is the fermion partial- wave operator. This latter operator 
should be constructed with the use of the orbital-momentum operator X*^^) and spin operator 
for fermion-antifermion system. 

We have two fermion spin states, 5" = and 5* = 1. For the spin-0 state, J = L, and for 
the spin-1 state, one has J = L — 1, L, L + 1. 

For fermion operator Q^^^l^^'^}^j_^^j-, one implies the same constraints as for the boson one: 
the fermion operator should be symmetrical, P-orthogonal and traceless: 

n(5,L,J) = n('5'^-^) p n(^'^.-^) =0 a C^'^'^"^) = f70) 

V/xi/x2.../ij_iAJj v^2Mi---M./-iMj ' MV^^2--./ij-iAtj ^' y^J.^^J.k^ ^J,l. ..fij \"^) 

The spin-0 operator for proton-antiproton system, F^^-*, is proportional to the 75-matrix. We 
normalize F*-°^ by the condition: 

Sp (F(0)(m + gi)F(o)(m - q^)) = 1 , (71) 

that gives 

p(o) _ _ .-^2) 

The angular-momentum operator for the spin-0 state is a product of F^°) and the angular- 
momentum operator Xjjp ^^{q): 

«id..-.M.(^) = r(°^4^,.(?)- (73) 

The spin-1 operator is constructed as follows: 

r(i) = ^ 1^7^ ^-^ ] =^(1^- ^" ] (74) 
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Here we have used that q = {ki — k2)/2 = m. The operator T^^^ is orthogonal to the total 
momentum P [PaT^^^ = 0) and normalized as follows: 

Sp (ri^) (m + q,)r^P (m - q^)) = gip . (75) 

The operators F^^^ and r*^"^ are orthogonal to one another in the spin space: it means that, if 
initial fermions are not polarized, the spin-0 and spin-1 states do not interfere with each other 
in the differential cross section. 

The spin-1 state with J = L — 1 is constructed from the operator F^^^ and angular- 
momentum operator X^;'^]} ,, . It carries J indices, so two indices, one from T^^^ and another 
from Xj/^^^Jj^j_^_^ , are to be absorbed, that can be done with the help of the metric tensor gafij+i '■ 

Qi\'^.i="-^n?) = r«4?..,,_,a(?) • (76) 

For J = L, the construction of the operator Q^}'^'^^^^ is performed by using antisymmetrical 



tensor e^jyii/ji/g: the operator F^^X^-^^^ must have the same number of indices as the angular- 
momentum operator. The only possible non-zero combination of antisymmetrical tensor, the 
X^"^^ and F*^^) operators, is given by the following convolution: 



e 



p (a) (77) 



However, the operator entering equation (77), being P-orthogonal and traceless, is not symmet- 



rical. The symmetrization can be performed by using the tensor Zj/^^ ^^^^ „ instead of Xj/J^ 
In this way, we get: 

^(l,L,J=L)/ N p pC^vC"') (^Q) 

Following the same procedure, we can easily construct the operator for the total angular mo- 
mentum J = L + 1. One has: 

Q£^.i="^^H?) = r^Z^^^Jq) . (79) 



After putting L = 1 in (|76|), (|78D and (ffOD, we have the P-wave operators which are used in 
the analysis of the pp and np annihilations. 

The moment-operator expansion presented above was used in analysis of the meson spectra 
in a number of papers |^, |ll|. 

Sometimes this technique is misleadingly referred as the Zemach expansion method. Com- 
paring the operators of Eqs. (|^) (or (|6ll)) and (H), one can see the common and different 
features of the three-dimensional approach of Zemach |^ and covariant method applied here 
and in p|, ||, |Tl|. For the operator (|6T1) used in the centre-of-mass frame, the expressions used 
in both approaches coincide. Indeed, the four-momentum has space-like components only, 
q-^ = (0,g), so the operator X^^^ ^ti(?) turns into Zemach's operator. However, for the ampli- 
tude (|65|) a simultaneous equality to zero of operators with zero components is impossible. In 
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59[| a special procedure was suggested for such cases, namely, the operator is treated in its own 



centre-of-mass frame, with subsequent Lorentz boost to a needed frame. But in the procedure 
developed in [§, |l^ and summarized in [^], these additional manipulations are unnecessary. 

The Lorentz boost should be also carried out upon the three-particle production amplitude 
considered in terms of spherical wave functions as well as in its version suggested by K 



Appendix C: Quark-combinatorics relation for the decay 
couplings 

In Table 7, we give the decay constants for the transition 

qq = nn cos (f + ss sin (f —* two pseudoscalar mesons , (80) 

where nn = [uu + dd) / -\/2, while the angle G defines the quark content of 77 and 77' mesons 
assuming them pure qq states: 77 = rzncosG — sssinB and 77' = nnsmQ + sscosG, with 
cos 9 ~ 0.8 and sin 9 ~ 0.6. 

The relations for the decay constants are given for planar diagrams which are the leading 
ones, in terms of the 1/N expansion ||3^, see also P, ^ for details. 



The gluball decay is the two-stage process: gluehall {qq)giuebau mesons. The transition 
gluehall — > {qq)giuebaii is not suppressed in the framework of expansion rule, see 0, p7[ , 
hence the transition constants 

gluehall two pseudoscalar mesons (81) 

are of the same order as the decay constants of the gg-state. The constants for the transition 
(pip obey the same relations as Eq. (75), with a special fixation of ip. By substituting 

^ ~^ '^gluehall , 90 ^ Gq , (82) 



we have tanipgiuebaii = y A/2. Such a definition of the (fgiuebaii is due to the equality {qq)giuebaii = 
{uu + dd + ^/Xss) I -\/2 + A: recall that we assume the new gg-pairs to be produced in the 
proportion uu : dd : ss = \ : \ : \. 

The coupling constants given in Table 7 satisfy the sum rule: 

gl{nn^c)h+ gliss^e)h = ^g',i2 + \). (83) 

c=inT,K K ,riri,T]r]' ,ri'ri' c=K K ,r]r],r]ri' ,ri'r)' 

The factor (2+ A) is due to the production of one additional gg-pair in the decay of the gg-meson; 
Ic is the identity factor, see Table 7. 

For the glueball decay the sum of couplings squared over all channels is proportional to the 
probability (2 + A)^ to produce two gg pairs. So 

E Gl{e)I{e) = \Gl{2 + \f. (84) 

c=iviT ,K K ,rir],T]r]' ,r]'rj' 
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Ill 



The relations for the decay couphng constants in case of non-planar diagrams may be found 
The analysis proved that the non-planar diagram contribution is suppressed — 



11 



just as it should be within the rules of 1/A^ expansion, so we do not use such type of terms in 
this analysis. 
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Table 7: 

Coupling constants given by quark combinatorics for a g'^-meson decaying into a pair of pseu- 

doscalar mesons in the leading terms of the 1/Nf. expansion. 99 is the mixing angle for nn 
and ss states, and is the mixing angle for rj — rj' mesons: rj = nn cos — ss sin and 
T}' = nnsin0 + sscos0. Glueball decay couplings in the leading terms of 1/iV expansion are 
obtained by the replacements gocoscp — > V^Gq, 5^0 sin <^ — > VXGo. 



Channel 


The gg-meson decay 

couplings in the 
leading terms oll/N 
expansion 


Identity 
factor in 
phase space 






go cos If/ V2 


1/2 






go cos (fi/ V2 


1 


K+K- 




go{V2sm(p + y/Xcos(p)/^/8 


1 






go{\/2sm(p + -s/A cos (^) /-s/8 


1 


VV 


9o ( 


cos^O cos</?/-\/2 + -s/A sin(/7 sin^O^ 


1/2 


rjrj' 




sin cos ^cos </?/ V2 — \/\ sin ip^ 


1 


rj'rj' 


9o ( 


sin^ cos(^/\/2 + \/A svnip cos^0j 


1/2 
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Figure 1: A mass projection of the acceptance-corrected Dalitz plot for the pp annihilation into 
7r°7r°7r° in hquid H2. The curve corresponds to the fit in Solution II-2. 
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Figure 2: Mass projections of the acceptance-corrected Dalitz plot for the pp annihilation into 
n^rjr] in liquid H2. Curves correspond to the fit in Solution II-2. 
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Figure 3: Mass projections of the acceptance-corrected Dalitz plot for the pp annihilation into 
TT^TT^r] in liquid H2. Curves correspond to the fit in Solution II-2. 
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Figure 4: a,b) Mass projections of the acceptance-corrected Dalitz plot for the pp annihilation 
into 7r+7r''7r~ in liquid H2, c) the angle distribution between charged and neutral pions in c.m.s. 
of 7r"'"7r~ system taken at masses between 0.95 and 1.05 GeV, d) the angle distribution between 
charged pions in c.m.s. of tt^tt" system taken at masses between 1.35 and 1.55 GeV. Figure 4d 
shows the event distribution along the band with the production of p(1450). 
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Figure 5: a,b) Mass projections of the acceptance-corrected Dalitz plot for the pp annihilation 
into 71^71^71^ in liquid D2, c) the angle distribution between charged and neutral pions in c.m.s. 
of tt^tt" system taken at masses between 1.20 and 1.40 GeV, d) the angle distribution between 
neutral pions in c.m.s. of 7r°7r~ system taken at masses between 1.35 and 1.55 GeV 
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Figure 6: a,b) Mass projections of the acceptance-corrected Dalitz plot for the pp annihilation 
into 7r~Tr~'K'^ in liquid D2, c) the angle distribution between charged and neutral pions in c.m.s. 
of 7r~7r~ system taken at masses between 1.20 and 1.40 GeV, d) the angle distribution between 
charged pions in c.m.s. of 7r~7r° system taken at masses between 1.35 and 1.55 GeV. 
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Figure 7: a,b) Mass projections of the acceptance-corrected Dalitz plot for the pp annihilation 
into KsKsTT^ in liquid H2, c) the angle distribution between kaons in c.m.s. of Kstt^ system 
taken at masses between 1.20 and 1.40 GeV, d) an angle distribution between kaon pion in 
c.m.s. of KsKs system taken at masses between 1.25 and 1.45 GeV. 
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Figure 8: a,b) Mass projections of the acceptance-corrected Dalitz plot for the pp annihilation 
into K'^K'T^^ in hquid H2, c) the angle distribution between kaons in c.m.s. of Kt^^ system 
taken at masses between 1.20 and 1.40 GeV, d) the angle distribution between kaon pion in 
c.m.s. of K^K~ system taken at masses between 1.25 and 1.45 GeV. 
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Figure 9: a,b,c) Mass projections of the acceptance-corrected Dalitz plot for the pp annihilation 
into KlK^ti^ {KlK^tt~) in liquid H2, d) KsK mass projection of the acceptance corrected 
Dalitz plot for the pp annihilation into KsK~'n:^ . This reaction has some problems with 
acceptance correction and was not used in the analysis. The full curve corresponds to the fit 
of pp — > KlK~t:^ reaction normalized to the number of KsK'~7r'^ events. 
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Figure 10: a,b) Mass projections of the acceptance-corrected Dalitz plot for the pp annihilation 
into KsKsTT" in liquid D2, c) the angle distribution between kaons in c.m.s. of KsTr~ system 
taken at masses between 1.20 and 1.40 GeV, d) the angle distribution between kaon and pion 
in c.m.s. of KgKs system taken at masses between 1.25 and 1.45 GeV. 
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Figure 11: a,b,c) Mass projections of the acceptance-corrected Dalitz plot for the pp annihilation 
into KsK~7r^ in liquid D2, d) the angle distribution between Ks and 7r° in c.m.s. of KgK" 
system taken at masses between 1.25 and 1.45 GeV. 
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Figure 12: Comparison of the GAMS and E852 data in the |t|-interval 0.3 < |t| < 0.4 GeV^. 
Full circles are E852 data and open circles correspond to the substruction of the two sets of 
GAMS data: A^[0.3 < \t\ < 1.0 GeV^]/{20 MeV) - N[OA < \t\ < 1.0 GeV2]/(20 MeV). 
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Figure 13: Description of the E852 data at different f-intervals. Dashed curve shows the 
contribution from Oi-trajectory and dotted curve from TTdaughter-^'^^i^ctory. 
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Figure 14: Description of the GAMS data at different t-intervals. The dashed curves refer to 
the solution in the previous analysis [11]. 
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Figure 15: The ^-dependence of the K-matrix couphngs. a) 7r-exchanges: full curve is for 
j6are(720) and dashed one for /o^"''^(1250), b) 7r-exchanges: full curve is for /o^«^^(1230) and 
dashed curve for /q"^^(1600). c,d) Oi-exchanges: ^-dependence for the same states as in figures 
a and b. 
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Figure 16: Complex M-plane: trajectories of the poles for /o(980), /o(1300), /o(1500), /o(1750), 
/o(1200 — 1600) during gradual onset of the decay processes. 
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Figure 17: Correlation curves on the {(f,G/g) and {ip^g/G) plots for the description of the 
decay couplings of resonances (Table 5) in terms of quark-combinatorics relations (22). a,c,e) 
Correlation curves for the q'^-originated resogig-nces: the curves with appropriate A's cover 
strips on the {ip,G/g) plane. b,d,f) Correlation curves for the glueball descendant: the curves 
at appropriate A's form a cross on the {(p, g/G) plane with the center near (/? ~ 30°, g/G 0. 




Figure 18: Linear trajectoties on the (n, M^)-plot for scalar resonances (a) and bare scalar 
states (b). Open points stand for the predicted states. 
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